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Abstract 

In this paper we consider the direct scattering problem of obliquely incident time-harmonic 
electromagnetic plane waves by an infinitely long dielectric cylinder. We assume that the 
cylinder and the outer medium are homogeneous and isotropic. From the symmetry of the 
problem. Maxwell’s equations are reduced to a system of two 2D Helmholtz equations in the 
cylinder and two Helmholtz equations in the exterior domain coupled on the boundary. We 
prove uniqueness and existence of this differential system by formulating an equivalent system 
of integral equations using the direct method. We transform this system into a Fredholm type 
system of boundary integral equations taking advantage of Maue’s formula for hypersingular 
operators. Applying a collocation method we derive an efficient numerical scheme and provide 
accurate numerical results using as test cases transmission problems corresponding to analytic 
fields derived from fundamental solutions. 

Keywords direct electromagnetic scattering, oblique incidence, integral equation method, 
hypersingular operator 


1 Introduction 

An interesting area of electromagnetism for its applications and the arising theoretical problems 
is the scattering process from obliquely incident time-harmonic plane waves by an infinitely long 
cylinder. The basic waves in the propagation domain satisfy Maxwell’s equations [f, 3, 16, 18] and 
due to the symmetry of the problem it is equivalent to find two scalar fields satisfying a pair of two- 
dimensional Helmholtz equations with different wavenumbers. The complication appears in the 
boundary conditions. Even for the case of a perfect conductor, in the boundary conditions appear 
tangential derivatives which make the analysis more difficult. There are many studies providing 
analytical or numerical solutions [2, 13, 19, 20, 21, 22, 24, 25]. The proposed methods are based 
on specific geometries or well known numerical schemes without examining the well-posedness of 
the corresponding boundary value problem. 

Recently, Wang and Nakamura [23] used a more elegant theoretical analysis to prove well- 
posedness of the problem based on the integral equation approach. They proved theoretical and 
numerical results for the case of homogeneous impedance cylinder using integral equations. For 
the theoretical analysis they used properties of the Cauchy singular integrals and proved that the 
derived system is of Fredholm type with index zero. For the numerical results they applied a specific 
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decomposition of the kernels and formulations using Hilbert’s and Symm’s integral operators. 
Considering trigonometric interpolation, they introduced an efficient numerical scheme. 

The case for general dielectric cylinders is not considered yet, however the same authors, in a 
later work [17], investigated a more complicated model having also a non-homogeneous part, in 
the sense that the permittivity and the permeability of the exterior medium are non-constants and 
smooth in a bounded domain surrounding the cylinder. The main theoretical analysis providing 
uniqueness and existence in non-homogeneous materials is much harder. For the well-posedness 
they used the Lax-Phillips method [7] . 

In this work, we examine the case of infinite dielectric cylinder illuminated by a transverse 
magnetic polarized electromagnetic plane wave, known as oblique incidence. More precisely, in 
the second Section starting from Maxwell’s equations we describe initially the derivation of the 
mathematical model for the scattering process from obliquely incident time-harmonic plane waves 
for the case of infinite inhomogeneous cylinder. We assume that transmission conditions hold on 
the boundary. The boundary conditions involve normal and tangential derivatives of the fields. 

In Section 3, we formulate the direct problem in differential form. We derive the Helmholtz 
equations and the exact form of the boundary conditions in the case of homogeneous cylinder. We 
prove that the problem is uniquely solvable using Green’s formulas and Rellich’s lemma. Consid¬ 
ering the direct method, initially applied in transmission problems in [5, 8, 9], we formulate the 
problem into an equivalent system of integral equations. We show that this system is of Fredholm 
type in an appropriate Sobolev space setting. Due to uniqueness of the boundary value problem 
existence follows from Fredholm alternative. The system consists of compact, singular and hyper¬ 
singular operators. We consider Maue’s formula [14], as in the case of the normal derivative of the 
double layer potential, to reduce the hypersingularity of the tangential derivative of the double 
layer potential. 

In the last Section we investigate numerically the problem by a collocation method based on 
Kress’s method for two dimensional integral equation with strongly singular operators [10]. We 
transform the system of integral equations to a linear system by parametrizing the operators 
and considering well-known quadrature rules. We derive accurate numerical results for the four 
fields, interior and exterior, and we compute numerically the far-held patterns of the two exterior 
helds computed for a specihc boundary value problem. Namely, we consider boundary data corre¬ 
sponding to analytic helds derived from point sources, where the interior and exterior helds have 
singularities outside of their domain of consideration. 


2 Formulation of the direct scattering problem for an inho¬ 
mogeneous cylinder 


We consider the scattering problem of an electromagnetic wave by a penetrable cylinder in R^. 
Let X = {x,y,z) S R^, then we model the cylinder as = {x : {x,y) G D, z € R}, where ft is 

a bounded domain in R^ with smooth boundary F. The cylinder flint is oriented parallel to the 
z-axis and fl is its horizontal cross section. We assume constant permittivity Cq and permeability 
y.Q for the exterior domain flext •= \ ^int- The interior domain flint is characterized by the 

electric constants y,{x) = n{x,y) and e{x) = e{x,y) for all {x,y) G D, z G R. 

We dehne for x G flext, t G R the magnetic held and electric held E^^*{x,t) and 

equivalently the interior helds H™^{x,t) and for x G flint, t G R. Then, these helds 

satisfy the Maxwell’s equations 


vx + ^o 
vx -k y 


dH^^xt 

dt 

dH^nt 

dt 



Vx - Co 

Vx - e 


d£;ext 

dt 

dE’-^t 

dt 



X G fl^xt, 
X G flint- 
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On the boundary F, we consider transmission conditions 

n X = h X h X = n X a; e F, 

where h is the outward normal vector, directed into riext- 

In order to take advantage of the symmetry of the specific medium, we probe the cylinder with 
an incident transverse magnetic (TM) polarized electromagnetic plane wave, the so-called oblique 
incidence in the literature. An arbitrary time-harmonic incident electromagnetic plane wave has 
the form: 

E^^%x,t;d,p) = T^T^Vx Vx 6"“*, 

V / 

W'^%x,t;d,p) = —-Vx 

where a; > 0 is the frequency, ko = is the wave number in the exterior domain, p is the 

polarization vector and d the vector describing the incident direction, satisfying d E p. 

In the following, due to the linearity of the problem we suppress the time-dependence and we 
consider the fields only as functions of the space variable x. In order to describe the incident fields 
for the specific TM polarization, we define by 0 the incident angle with respect to the negative x axis 
and by (f) the polar angle of d (in spherical coordinates), then, d = (sin0cos<(), sin0sin<(), — cos0) 
and p = (cos0cos()), cosSsini^, sin0), assuming that 0 G (0,7r/2) U (7r/2,7r). Hence, we obtain 

E^^^(x; d,p) = ^dxpx de^’^°^-^ = ^ 

ff“^(a;;d,p) = -^d x ^(sin<(), - cos 0) 

Vih) 

Taking into account the cylindrical symmetry of the medium and the z-independence of the 
electric coefficients we express the incident fields as separable functions of {x,y) and z. Thus, we 
define P = ko cos 9 and kq = ^/k^ — = ko sin 9 and it follows that the incident fields can be 

decomposed to 

£;“^(a;;d,p) = e“=(x, y) W'^%x-d,p) = y) (2) 

where 

e™'^{x,y) = ^ ^(;»«o(a:cos0+i/sin<ii) 

Ii“"(a;,y) = ^ (sin (/), - cos (^, 0) 

V/^ 

Now, we are in position to transform equations (1) into a system of equations only for the 
z-component of the electric and magnetic fields. Firstly, we see that for the specific illumination 
of the form (2), using separation of variables, also the scattered fields take the form: 

£;®^(a;;d,p) = e'*^(a:,j/)e“*^^, iT®^(a;; <i,p) = ft,^°(a;, y) x G ^lext, 

where = (ef^, e|'^, e§'^) and hf^). Then, the exterior fields are given by 

E--\x;d,p) = {e^%x,y) + E^%x,y)) = e^^\x,y) x G 

H^-*(x;d,p) = (h^^(x,y) + h^^^(x,y)) = e«"‘(cr, y) x G 
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Equivalently, the interior fields are represented by 

E^^\x-d,p) = W^\x;d,p) = h^^\x,y) x e 

where e“* = (ef*, e|”*, and 
For any field of the form 

E{x\d,p) = e{x,y) , H{x;d,p) = h{x,y)e~^^^, a; e R^, 


we consider the Maxwell’s equations in R^ for arbitrary e, y and = ytoo^ — /3^ (remark here the 
space dependence of e,y). Then, following [17] we obtain the relations 


ei{x, 

y) = 

F 1 

- iyuj 

dh-i 

dy 

62(2;, 

y) = 

-fI 

+ iyuj 

dhs 

dx 

hi{x 

,y) = 

1 

f n '9^3 

des 

~¥ 

\ip — [x,y) +ieuj 

dy 

h 2 {x 

,y) = 

1 

( '5^3 

de^ 

~¥ 

hp-^(a:,y) - iew 

dx 


( 3 ) 


Substituting (3) in (1), we have that the pair (e^,h^) satisfies the equations 


— V-(^^ 
euj 


-V' ( —:^Vh-3 I — - J \I I — 1 • 

yuj \ 


where 


pui 


J = 




Ves + k'^h'i = 0, 


0 1 

-1 0 


The interior and the exterior domains are characterized by different wavenumbers, given by 


k‘^{x) = 


■= y) e[x, y) u? - 0 ^, xG 


klAx) := yoeoUJ^ - = 




X gVL^ 


In this section, for completeness in the formulation of the direct problem we keep the space 
dependence of kint- Later, we consider only the case of constant parameters. Here, we have to 
assume that p{x)e{x) > eoPo^osO in order to have inf,j. kf^f.{x) > 0. Thus, the fields el^*{x,y) and 
hl^*{x,y) satisfy 


ef * = 0 , A/if * + = 0 , xG^^^t, 


( 4 ) 


and the interior fields 


kfntjx) 

e{x) 

kLtjx) . 

fj.{x) 


e{x) 


Vel 


kfnti^) 

viuint 

,2 


+ ^^JV 


kfntix) 


e{x)oj 

yL{x)uj 




JV 


kfnti^) 

P 

kLti^) 


Whl-* + klPx)eT^=0, xGn, 


■Vef‘ + AL(a;)/ir‘=0, xGn.^t. 


( 5 ) 


Now, we are going to derive the exact form of the boundary conditions. We introduce the 
notations: = x ei + y 62 , ht = x hi + y / 12 , and V* = x^ + y-^, where x, y, z denote the unit 
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vectors in R^. Let (n, r, z) be a local coordinate system, where fi = (ni, 712 ) is the outward normal 
vector and t = (—n 2 ,ni) the outward tangent vector on F. Then, from (3) we obtain 


T • e* = -p- {ipLwfi ■ Vths + iPf ■ VtCa), 

T • /it = -p {-ieujh ■ VtCa + i/ 3 f • Vt/ 13 ), 

using that t ■ (z x Vf) = n • Vf • . 

We observe, setting zero to the z-component of n, f in R^, that 

n X E = —ear + (11162 — 11261 ) i, fix H = —h^f + (ni /12 — 112 / 11 ) z. 

Then from (3) and (6), we derive 


h X = -ef‘f + f • er* i, 


fi X + f • /i“*, 


Cg / “T » -^5 /fc A -t-t — /i'3 

for the exterior fields, where := ® 6“‘ + y e“‘, := & /i“* + y and 


( 6 ) 


fi X = -e“*f + f • el 


fi X + f • fi™* i. 


for the interior helds, where e)”‘ := + ye™*, /i™* := x/i™* + yh™*. 

Here, we observe that the tangential forms of the fields can be written in terms of f and z, two 
linear independent vectors. Thus, the boundary condition 


fi X £;“* = fi X a;er. 


is equivalent to the system 


-3 — 

and equivalently for the magnetic fields 


e™* = 6^^‘, f • e™‘= f • e“‘, a; e T, 


We define 


= f-fi™* = f-/i®^‘, a;er. 


9 . ^ 9 ^ 

— =n-Vt, ^ = "T • Vt 

on OT 


and we rewrite the above boundary conditions as 


int _ ^Gxt 


and 


to = e 


3 ^ 


9n ' kf^tix) 9 t dn 


y{x) 9 / 13 ”* /? 963 ”* ^0 9hl 


-UJ- 


+ 


+ 


■Lint _ j.^ext 


dr 


t{x) 13 cq /3 dhf^ 

' ' — — 


a: e r, 
a: G r, 

a; G r, 
a; G r. 


(7) 


( 8 ) 


9n fcL(a=) 9t dn dr 

To ensure that the scattered fields are outgoing, the components must satisfy in addition the 
radiation conditions in R^ : 


5e? 


dr 


f 9 hl 


lim yfr ( -) = 0, lim ( - 1 ^ -) = 0, 


V dr 


(9) 
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where r = |(a:, j/)|, uniformly over all directions. 

Thus, the direct transmission problem for oblique incident wave, is to find the fields h™*, e™* 

and e|'^ which satisfy equations (4) and (5), the transmission conditions (7) and (8) and the 
radiation conditions (9). 

We remark here that since we consider TM polarized wave, see equation (2), the incident helds 
for X € ^ext are simplified to 

ef^{x,y) = hl^^{x,y) = 0. (10) 


3 The direct problem for a homogeneous cylinder using the 
integral equation method 

From now on, x G R^. In this section we consider the simplified version where y,{x) = and 
e{x) = ei are constant in the interior domain. To simplify the following analysis, we set C 

R^, flo = R^ \ fl and 

uo{x) = ef{x), vo{x) = hf{x), x G ^o, 
ui{x) = e“*(a;), ui(a;) = h™*{x), x G fli. 

In the following, j = 0,1 counts for the exterior {x G fig) and interior domain (x G fii), re¬ 
spectively. Then, the direct scattering problem, presented in the previous section, is modified to 


Auj + Kj Uj = 0 , Avj + Vj = 0 , a; € rij, 
for j = 0,1 where nf = yieiui'^ — j3‘^, with boundary conditions 

a: e r, 

a: e r, 
a: e r, 

a: G r, 

where jlj = yj/Kj, ej = Pj = P/kj and the radiation conditions 


ui=uo + eg"'", 

. dvi . o 9ui _ ^ , p ^'“0 , « 

-^ - f Po-^ - f — I 

on OT on OT or 


Vl = Vg, 


. dui dvi 

- 

on OT 


_ - dug dvg 

~ -^ -Po^^i 

on on OT 


lim ^/r ( - IKgUg 
r—>-oo y or 


= 0 , 


lim a/t f ^ - iKoUo ) = 0. 

7 —yea \ OT J 


( 11 ) 

(12a) 

(12b) 

(12c) 

(12d) 

(13) 


Theorem 3.1 If k\ is not an interior Dirichlet eigenvalue, then the problem (11) - (13) has at 
most one solution. 


Proof: It is sufficient to show that if ug,vg,ui,vi solve the homogeneous problem (11) - (13), that 
is for Cg”'^ = 0, then, uq = I'd = 0 ia and ui = Ui = 0 in fli. Let Sr be a disk with radius r, 
boundary T^, centered at the origin and containing fli. We set Itr = Sr\ fli, see Figure 1. 

The boundary conditions of the homogeneous problem read 


Ml = Mo, 


_ dvi _ dvg Pi dui 

Pi o Po o o 

on on UJ OT 


Pg dug 

UJ dT ’ 


Vl = Vg, 

_ 9 mi _ dug Pi dvi Pg dvg 
Cl eg „ ... .. , 

On on UJ OT UJ Ot 


X GT, 
X GT, 
X GT, 
a; G F. 


( 14 ) 
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Oq 



We apply Green’s first identity in and considering (11) we obtain 


ei [ Ui^^ds = ei [ (iVuiP + MiAmi) da: 
Jr dn 


= h j (|Vuip - Kiluip) dx, 
Ml / vi'^ds = jli / {iVvil"^+ viAvi) dx 


dvi 


= fj-i [ (iVmp - Kilmp) dai. 

J Oi 

Similarly, Green’s first identity in together with equations (14) and (15) gi 

eo J uo^^ds = eo J {\Vuo\‘^ + uqAuo) dx + cq J ds 

= eo (|VuoP - KqImoP) rfa: 
jQr- 


gives 


. dui Pi dvi Po dvo 1 , 

“0 ei ^-^- 1 -q— I ds 

on uj or Lu or 


= eo (|VuoP - KqImoP) da; 
jQr- 

f /|^ ,2 2, |2\ , /3l f a , /3o /■ dvo 

dni ' UJ Jr dr uj Jr dr 


(15) 
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and 


dvp 

dn 


dvp 

dn 


A*o / '^o^^ds = ij.o J (|Vvo|^ + ?^oAvo) da: +/io j vo^ds 

= Mo/ (|Vwo|^ - KokoH da; 

Jn^ 

dvi Po dup Pi dui 


+ / ^'o Ml 


dn UJ dr UJ dr 


= pp (|Vz;oP - KokoP) da; 

+ Pi [ (iVnip - Kiliiip) da; -— 

JQo ^ 


ds 


dup 

vp ds 


We add the above two equations and noting that 


Pi 

UJ 


dui 

vi^ds. 


dv 


dui 


dvp 


dup 


— / ui——ds = / —ds, / up-^ds = — Vp— —ds, 


we obtain 


dr 


3 m I eo 


dr 


dup 

Up ds + fxp 
dn 


Vp 


dr 
dvp 


dr 


dn 


ds ) =0, 


or equivalently, using the radiation conditions (see equation 2.12 [23]) 


lim 

r—>-oo 



^eoluoP + Co 


dup 

dn 


2 

+ MokoP + Mo 


dvp 

dn 


^ ds = 0. 


Thus 



lim / jzJopds = 0, 
r^ooJr^ 


and by Rellich’s lemma it follows that mq = Uq = 0 in fig. Hence, rtg = Ug = 0 in T and ui = Ui = 0 
in r from the boundary conditions. Then, rti = ui = 0 in Hi follows from the unique solvability 
of the interior Dirichlet problem, given the assumption of the theorem. □ 


We define the fundamental solution of the Helmholtz equation in : 

^j{x,y) = - y\), a;,yeHj, x^y, (16) 

where idg^^ is the Hankel function of the first kind and zero order. For a continuous density /, we 
introduce the single- and double-layer potentials defined by 


iSjf){x) = J $j(a;,y)/(y)ds(y), x G H^-, 

f f)^ ■ 

iDjf)ix) = ^^^(a;,y)/(y)ds(y), x G %, 


(17) 
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and their derivatives, normal and tangential as a; —>■ F, using the standard jump relations, see for 
example [3, 6] 


d r (9$, 1 1 

^('S'j/)(3:) = i^^^{x,y)fiy)ds{y) t ^/(x) ■= (NSjf){x) t ^^(x), ® G T, 

d f 

^{D,f){x) = Q^^J±g^i^^y)ny)ds{y) := {ND,f)ix), x G F, 


i^SjfKx) = j^^^{x,y)f{y)dsiy) 


iTS,f){x), 


£c e F, 


dr 


{Djf){x) 


dr(t)Z(y) ^^’ yVMds(y) ± i (X) := (TD,f)M ±\%{<c). x e T, 


( 18 ) 

where the upper (lower) sign indicates the limits obtained by approaching the boundary F from 
rio (fli), this means when j = 0 (j = 1). For the last equation, we assume / to be continuously 
differentiable. We presented the jump relations for continuous densities only for simplicity, since 
in the following proof we use a Sobolev space setting. 

Flere, we have to mention the continuity of the single-layer potentials in and the discon¬ 
tinuity of the double-layer potentials (±1/). All the integrals are well defined and particularly 
the potentials Sj,Dj and NSj have weakly singular kernels (logarithmic singularity), the poten¬ 
tials TSj have Cauchy type singularity (of order l/\x — y\) and the potentials NDj.TDj are 
hypersingular (of order l/|a: — yp). 


Theorem 3.2 If is not an interior Dirichlet eigenvalue and Kq is not an interior Dirichlet and 
Neumann eigenvalue, then the problem (11) - (13) has a unique solution. 


Proof: We apply the direct method, see for instance [9], to transform the problem into a system 
of integral equations. We consider the Green’s second theorem in the interior domain 

-ui{x) = [x,y)ui{y)ds{y) - J^^i{x,y)J^^{y)ds{y), (19) 

= {Diui){x) - {Sid^ui){x), xGfli, 


similarly 



-vi{x) = {Divi){x) - {SidrjVi){x), 

X G rii. 


and in the exterior domain 



uo{x) = {DoUo){x) - {Sodr,uo){x), 

X G flo, 

(20) 

vo(x) = (DoVo)(x) - (Sod,jVo)(x), 

X G fig- 

Letting a; — )■ F, in the above formulas and taking the normal and the tangential derivatives 

on F, 


we obtain 


{NS, ± i/) d,u, = ND,t 
{ ddj T Uj = Sjd, 




TDjUj — TSjdrtUj = -p^drUj, 


{NS,±\l)dr^Vj = NDjV,, (21a) 

{dd^i \d) Vj = SjdjjVj, (21b) 

TDjV, - TSjdnVj = ±\drV,. (21c) 


Combining the relations in (21b) for j = 0 with the boundary conditions (12d) and (12b) respec¬ 
tively, we have 



MO = - Sod^er, 

eo eouj eouj 

Vo = ^Sod^vi + -^Sodrui - -i^SodrUo - -i^Sodref^ 
ho Mow /Tow /ioW 


( 22 ) 
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We define 


K,:={NS,±^I) ^ND„ L, := 2{TD, - TS,Kj). (23) 

The operator Kq is well-defined since the integral equation (21a) for j = 0 corresponds to the solu¬ 
tion of the interior Neumann problem and we assumed Kq not to be interior eigenvalue. Similarly, 
Ki is well-defined if ref is not an interior Dirichlet eigenvalue. 

Then, the system of equations (22), using (21), is transformed to 

(do- l-l) uo - - ^SoL^vi - ^SoLqvo = Sod^er, 

V 2 / eo eouj 

(dq — -l\ vq — ^SqKiVi + -^—SqLiUi + 

V 2 / fio Mow Mow Mow 

We consider now equations (12a) and (12c), to obtain 

(do- b) uo - "^SoKiuo - ^SoLivo - ^SoLqVo = Sodr^er + "^SoK^er, 

V 2 y Co CqW CqW eo 

(do- ll) VO - ^SoK^vo + ^SoL^uo + ^SoLqUo = -^Sod^e^ - ^SoL^e^- 
\ 2 / Mo Mow Mow Mow mow 

The above system in compact form reads 

{D + K)u = b, (24) 


where 


D = 

K = 


u = 


0 Do-\lj 

-hSoK_ 


-j^SoiPiL^+/3oLo)\ 


^S'oitTi 

fJ-0 ^ ^ 


+l3oLo) 

( -Sodr, + %SoK, 

\-fk^So{Podr+PlL,) 




inc 
^3 • 


We assume that T is of class (7^’“, 0 < a < 1. We know that Do : H ^/^(T) — H ^/^(T) is 
compact, thus {Do — \l)~^ '■ —>• is bounded, if Kq is not an interior Dirichlet 

eigenvalue. Then, D : (7M“^/^(r))^ —>■ (7J“^/^(r))^ is bounded and (24) is transformed to 


{I + D-^K)u = D-^b. 


(25) 


First, we show that K is compact. We recall that 

SoKi = So{NS^-\l)~^ ND^, 

SoiPiL^ + PoLo) = 2So [PiTDi + PoTDo - PoTSo {NSo + i/)”' NDo) 
- 2piSoTSi {NSi -\l)~^NDi, 


and the following properties, see [4, 5], So '■ H ^/^(T) —H ^/^(T) is compact, NDj, TDj : 
iJi/2(r) ^ £f-i/2(r) are bounded, TSj : H-^/'^{T) £f-i/2(r) are bounded, and {NSj ± \l)~^ : 

^-i/2(r) aj-g bounded due to compactness of NSj : —>■ Then, 

the operators 

SoK^, So{PiL^+PoLo) ■■ ^ H-^/^{V), 
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are compact. Hence, K : {H is also compact resulting to the compactness 

of D-^K : (iJi/2(r))2 ^ (H-i/2(r))2. 

Next we prove the uniqueness of solutions of equation (25). Solvability follows from the Fred¬ 
holm alternative theorem. Let be the solution of the homogeneous form of (25). Then, 

the potentials 

uofi{x) = {Do(j)o){x) - {Sod^(j)o){x), x e Hq, 
vo,o(x) = {DoiIjo){x) - {Sodr,tpo){x), X e Ho, 

and 

Uia{x) = -{Di(l)i){x) + {Sidn4'i){x), x S Hi, 
ui,i(a;) = -{Ditlji){x) + {Sidr,tpi)ix), x G Hi, 

solve the homogeneous form of the problem (11) - (13). From Theorem 3.1, we have that uo,o = 
Wo,o = 0, in Hq and ui i = Ui,i = 0, in Hi where the densities depend on the solution of the 

homogeneous case. 

We construct 

uo,i(a:) = {Do(l)o){x) - {Sodrj(l)o){x), x G Hi, 
vo,i{x) = {Do'tljo){x) - {Sodr,tpo){x), X G Hi, 

and 

uifi{x) = -{Di(l)i){x) + {Sidn4'i){x), x G Hq, 
ui,o(a:) = -{DiiIji){x) + {Sidr,tpi){x), x G Hq. 

Considering the jump relations, at the boundary F we obtain 

“0,0 — Uo,l = <(>0, Ui^i — Ui^o = (j)i, 

“ 0,0 - “ 0,1 = ^/’ O , “ 1,1 - “ 1,0 = ' 01 - 

Since Uo.o = “i,i = “o,o = “i.i = Oj 't’o = 4>ii ipQ = V’l: on b, we find 

“0,1 =“1,0, “0,1 = “1,0, a; G F. 

Similarly, we can rewrite the other two boundary conditions of (14) for those fields taking the 
differences of the normal and tangential derivatives as a; —> F. Thus, we see that “o,i, “o,i, “i,o ^nd 
ui.o solve the homogeneous problem but with ki and kq interchanged and from Theorem 3.1 we 
get also uo,i = “o,i = 0, on F and hence (po = = 0 from (26). □ 

In order to handle the hypersingularity of the operators TDj we work in a similar way as 
Mitzner [15] derived the Maue’s formula [14] of the hypersingular operator NDj, namely 

{NDjf){x) = {x) + h{x) ■ {S, nf){x), xGT. (27) 

This transformation reduces the hypersingularity to singularity of Cauchy type (first term) and to 
a weak singularity (second term). 

Theorem 3.3 Let f G (^^’“(F). The hypersingular operator TDj can be transformed to 

iTD,f)ix) = -(^NS,^yx) + ^^fix)-iS,hf)ixf xGF. (28) 

Proof: We recall equation (19). Applying Green’s first theorem to v ■ Vui,i; arbitrary constant 
vector, and <i)i, using that 


A$i(a;,y) Ki$i(a;,y) = -(5(ai - y). 
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yields 


/ Vy{v-VyUiiy)) ■Vy<^i{x,y)dy - kI ^i{x,y) v ■ VyUi{y)dy 
J 


(29) 


f 

= v-\/^ui{x) + J ■^-^{x,y)v ■ \/yUi{y)ds{y). 

The first integral can be transformed to 

/ V(i) • Vui) • V$i dy = — / + V X (v X Vui)) • V$i dy — • / wiV^idy 

</ J Oi 

— V ■ y (n X V<i)i) X Viti ds{y). 

Then, (29) reads 

V ■ j Vy($i ui)dy + J {fix Vy$i) X V^ui ds(y) + J ds{y)^ = -v ■ V^ui. 

Using some vector identities and suppressing the inner products with v (holds for any vector), 
we end up to 


-V^ui = 


J X {fly X WyUi) + Kf^iuifiy + ds{y), 


for x G Qi. We multiple this equation with fi{x) (inner product) and considering the limit as x 
approaches the boundary T from inside and the corresponding jump relations, we obtain [15] 


l_d_ 
2 dn 


Ui = y ((fia; X Va;$i) • {fly X \7yUi) + Kl^iUi{fi^ ■ fly)) ds{y) 


f dui 
/r dn^ duy 


ds{y), £c G r. 


(30) 


Now, equating the above equation and the normal derivative of (19) as a; —)■ T we obtain (27) 
in R^. 

We take the tangential derivative of (19) and considering the jump relations we get 


d 


1 dui 


- — ui{x) = {TDiUi){x) - -^{x) - TSi— {x), a; e r. 


dr 


2 dr 


dui 


dn 


We replace fi{x) by f{x) in (30) (considering the appropriate jump relations) and restricting 
ourselves in R^, we have 

--—ui= J (-(rij, • Vx$i)(Ty • VyUi) + • fij,)) ds(y) - y ^^ds(y), a: G T. 

Observing the last two equations, we obtain (28), the equivalent of the Maue’s formula for the 
tangential derivative of the double-layer potential which also reduces the hypersingularity of the 
potential. □ 


4 Numerical results 

In this section we present numerical examples by implementing the proposed method. We use 
quadrature rules to integrate the singularities considering trigonometric interpolation. Regarding 
the convergence and the error analysis of the quadrature formulas, we refer the reader to [11] for 
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the weakly singular operators and to [10] for the hypersingular. We solve the system of integral 
equations considering these rules by the Nystrom method. 

We assume the following parametrization for the boundary 


r = {z{t) = (zi(t),Z 2 (<)) : t e [0,27r]}, 

where 2 : : R —> R^ is a C^-smooth, 27r-periodic and counter-clockwise oriented parametrization. 
We assume in addition that 2 : is injective in [0, 27r), that is z'{t) ^ 0, for all t G [0,27r]. 

Now, we transform the operators in (17) and their derivatives, see (18), into their parametric 
forms 

r 2 Tr 


Jo 

Jo 

p27T 

{NSji;)it)= M^^^{t,s)ilj{s)ds, 
Jo 

p2-n 

Jo 


and the special forms 




47r 


cot 


s —1\ dtp 


is)-M^^^it,s)tPis) 


ds 


r2TT 


+ k'^ {n{t) ■ h{s))M''^^{t,s)tp{s)ds, 

Jo 

{TDjtP){t) = (t, s)tP{s)ds 

p2-K 

+ {T{t) ■ n{s))M^^ {t,s)tp{s)ds, 

Jo 

for t G [0, 27r], j = 0,1, and tp{t) = f{z{t)), n{t) = n{z{t)), T{t) = f{z{t)), where 

= ^F^^^(Kj|r(t,s)|)|2;'(s)|, 


(t, s) = 


iKj n{s) ■ r{t,s) ( 1 ) 


4 |r(t,s)| 

iKj n{t) ■ r{t,s) ^( 1 ) 
^4 -|r(f,s)| 1 


i7r^(«:,|r(t,s)|)|^'(s)|. 


(t, s) = -^s)I)l^'(s)I, 






- 2k: 


H[^\Kj\r{t, s)|) 


iKi z'(t) ■ z'(s) , , , ,,, 1 . _9 


r{t,s)\ 
t — s 




- 2k. 


H[^\Kj\r{t,s)\) 

|r(t,s)| 


iK, Jz'{t)-z'is) ^,) 

+ 4 («,l^(i,s)|), 


(32a) 


(32b) 


(33) 
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with r{t, s) = z{t) — z{s), and 


{z'jt) ■r{t,s)){z'{s) ■r{t,s)) 

(J z'jt) ■ r{t,s)){z'{s) ■ r{t,s)) 

|r(t,s)|2 

Here, we have used the formulas (t) and (t) = The form 

in equation (32a) is based on (27), derived by Kress [10] and improved in [12], The derivation of 
(32b) is easier. Namely, we define dtj '■= J z'd^ and 

then 


M{t, s) = 
M'\t,s) = 


{z{t)) 


r27T 


Jz'{t) d , , z'is) d'if) . ,, ,, ,, , 


l^'WI 


1 ri r) 

w\i So" 

(t, s)tjjis)ds 


>0 

r2^ q2 


dtjds 


and (32b) follows by simply adding the parametrized form of the second term in the right-hand 
side of (28). The kernels in (33) admit the decomposition 


M'=(t,s) = 






for k = Sj,Dj,NSj,NDj,TDj where and are analytic, due to logarithmic singularity of 
the functions at t = s. The case of has to be treated differently because of the Cauchy type 

singularity of the kernel as t = s. Thus, we split the kernel as 

{t, s) = (t, s) In ^dsin^ ^ 
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The kernels are defined for i ^ s by, see [10] 




Mr^{t,s) 


-7zMi^j\r{t,s)\)\z\s)\, 


1 

47r 

Kj n{s) ■ r{t, s) 
A-k |r(f, s)| 


Mi^j\r{t,s)\)\z'{s)\, 


Kj n{t) ■ r{t, s) 
An |r(f, s)| 

Kj f{t) ■ r{t, s) 


An |r(f, s)| 

Kj z'{t) ■ z'{s) 
An |r(i, s)| 

1 

An 


Mi^j\rit,s)\)\z'{s)\, 

Ji{Kj\r{t,s)\)\z'{s)\, 


K'^jJo{Kj\r{t,s)\) - 2Kj 


Ji{Kj\r{t,s)\) 


Ji{Kj\r{t,s)\), 


— K‘^jJoiKj\r{t, s)\) — 2 k 


|r(t,s)| 

Ji{Kj\r{t,s)\) 


|r(f,s)| 


Kj J z'(t) ■ z'{s) 
An |r(f, s)| 


Ji{Kj\r{t,s)\), 


with diagonal terms 

= 0 , 

Mf {t, t) = 0 , {t, t) = 0 , 

= 0 , 

where Jq, Ji are the Bessel functions of order zero and one, respectively. The kernels , for 
t s, are given by 

M^{t,s) = M^{t,s) - Mf(t,s) In ^dsin^ (“^)) 

and 

{t, s) = {t, s) - {t, s) In ^4 sin^ (“^) ) “ ^ ’ 
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with diagonal terms 




i C 1 , /'ti', M 

1 h{t) ■ z''{t) 

47r |z'(t)| 

1 ri{t) ■ z"{t) 

47r \z'(t)\ 

1 nt)-^"{t) 

47r \z'{t)\ 


W(t)l 


M. 


ND, 


{t,t) = 


+ 




TTi — 1 — 2C — 21n (^'^\z'(t)\^ 


{z'{t)-z”{t)f z'{t)-z"'{t) 

47r|2;'(t)|'^ 127r|2'(t)p 

{z’{t)-z’'{t)){Jz'{^-z’'{t)) 
A'K\z'{t)\^ 


s^\z'{t)r 


1 


J z\t)-z'"{t) 
^ 127r|2'(t)|2 


where C is the Euler’s constant. For the last approximation, we used same arguments as in the 
case of . 

Considering the equidistant points tj = jV/n, j = 0, ...,2n — 1, we use the trapezoidal rule to 
approximate the operators with smooth kernel 


(.274 


'ip{s)ds 


2n-l 

IY. 

j=0 


and the following quadrature rules for the singular kernels 


^ In ^4sin^ Y 

1 ^27r / ( \ o 

4^ I (^) 

274 / a — f\ 2n-l 

cot f j V’(s)ds « 

^ ^ i=o 

with weights 

Rj^\t) = — — ^ — cos {m{t — tj)) —^ cos{n{t — tj)), 

m—1 

1 1 
= - ^ XI “ 4 cos(n(t - tj)), 

m=l 

^ ^ cos(nt)] cot ^ ^ 

Then, the system (24), or similarly (25), considering the above parametric forms of the integral 
operators and the quadrature rules, is transformed to a linear system by applying the Nystrom 
method. 

To illustrate the efficiency of our method, we consider two different cases. In the first example, 
motivated by [23], we construct a model where the scattered fields can be analytically computed 
and in the second one we consider the scattering of obliquely incident waves. 
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In both examples, the parametrization of the obstacle is given by 


z{t) = (2cost + 1.5cos2t — 1, 2.5sint), t G [0,27r]. 


4.1 Example with analytic solution 

We consider four arbitrary points 2 : 1 , 2:2 G and 2 : 3 , 2:4 G rig and we define the boundary functions 
fk, /c = 1,2,3,4 by 

/i = - H^o^\Ko\ri{x)\), 

r ~ „(!)/ I / ■ r^x) ( 1 ) f(a;) • r 3 (a;) 

/2 = -fJ-iUJKiHl '{Ki\r4{x)\) —^- I3iKiHI ^(ki r3(a;) )—^— 

\ri[x}\ r3(a;)| 

+ ^lQU}KoH[ \Ko\r2{x )\)+ ^oKoH[ \Ko\ri{x)\)-^ — - 

\r2[x)\ |n(a^)| 

/a = H^^\Ki\r4ix)\) - H^^\Ko\r 2 ix)\), 

r ~ „(!)/ I / )n’n(x) ■ r 3 (x) ( 1 ) f(a;) • r4(a;) 

f4 = -eiUJKiff^ '{m\r3{x)\) —hTTmi -4Ki|r4(a;)|)- 

+ eoWKoilj^^(Ko|ri(a;)|)- 


|r3(a:)| |r4(a;)| 

n{x)-ri{x) ( 1 ) f{x)-r 2 ix) 

■ - PoKqHI '(Ko\r2[x)\)- 


\ri{x)\ 


17-2 (a:) I 


where rk{x) = x — Zk- Then, the fields 

r(i) 


uoix) = H^"’{ko\x- zi\), Vo{x) = H^^\ko\x-Z 2\), XGfl 


Vi{x) = Hlf^\Ki\x-Z 4 \), ajGOi, 


Ui{x) = H^^\ki\x- Zsl), Vi{x) = 
solve the following problem 

Auq + Kg Ug = 0, Aug + Kg Ug = 0, X G fig, 
Awi + K^ui = 0 , Aui + K 4 ui = 0 , X G Qi, 


with boundary conditions 


Ui = Mg + /i, 


a; G r, 
/2, a; G r, 

7^1 = 7^0 + /a, a; G r, 

i9mi o 9vi _ ^ duo „ ^Mg 
--^ 0 “^- ^/4) 

or on 


. i9mi , „ i9ui . dvo , „ duo 
Hiu}— + Pi^— = + Po^— 

on OT on OT 


dn 


dr 


and the radiation conditions 


lim y/r ( - iKoUo ) = 0, 

r —>^00 \ OV 


lim y/r ( ^ - ZKgUg ) = 0. 


For this problem, we can derive again a system as (24), where now b is replaced by 

'-I^Sof4 + t^SoK4f4 + itSoLj3\ 

.-^Sof2 + >tSoKj3 - h-SoL^h ■ 


bf = 


(34) 


Given (34) and the asymptotic behaviour of the Hankel function [4] , we know that the far field 
patterns of ug and vq are given by 


^(x) = 


_4je"/4 

x/STTKg 


^(X) = 


-4ie^^/4 

V'Sttko 


X G S, 


(35) 
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Figure 2. The parametrization of the boundary P and the source points. 



.,00 

Uq 

,,00 

^0 


0 

7r/2 

TT 

0 

7r/2 

TT 

32 

0.006470 

0.009056 

0.003944 

0.004714 

0.007460 

0.001767 

64 

0.003157 

0.004729 

0.001820 

0.002249 

0.004288 

0.000601 

128 

0.001581 

0.002374 

0.000908 

0.001125 

0.002161 

0.000292 


Table 1. Absolute errors of the far field patterns of uq and vq for different orders n 
at discrete points t. 


where S is the unit ball. Numerically, the far field patterns are given by 

piizlA /*27r 

^ J [-iKoix ■ h(s))(po(s) - (Xo<^o)(s)] \z'is)\d 

piizlA /*27r 

z;o~(*) = . n{s))Ms) - {Koi’o){s)] \z'{s)\d 


s, 

S, 


(36) 


where (f := {(po, V'o)^ solves 

iD + K)<p = bf. 

Here, we have used the representations (20) for the exterior fields and the asymptotics of the 
Hankel function. The operator Kq is given by (23). 

We consider the points Zi = (0.5, 1) and 2:2 = (0, —0.5) in Hi and the points 2:3 = (1, 2) and 
2:4 = (0, —2.5) in Ho, see Figure 2. We set w = 1 and n = 32. The exact values (35) and the 
reconstructed (36) for (ei,/ri) = (3, 2) and (eoiMo) = (1; 1) are presented in Figure 3. The results 
are presented for 9 = 7 r/ 3 . In Table 1, we provide the absolute errors of the far field patterns for 
different values of n and t. 

As a general comment we could say that the reconstructions are accurate and illustrate the 
feasibility of the proposed method. However, the convergence is slower compared to the impedance 
cylinder case [23]. The main reason is the complexity of the matrix K involving the product of four 
operators: SqTSj {NSj + p) NDj,j = 0,1 resulting to an increase of the condition number. 
As 9 7 r/ 2 , the results improve considerably. 
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Figure 3. The far field patterns: Reconstructed (blue open circles) and exact (red 
solid line). 


4.2 Example with oblique incidence 

In this example, we consider the usual obliquely incident (TM) polarized electromagnetic plane 
wave, resulting to the forms (10). We keep the same values for all the parameters as in the previous 
example. We restrict the computations of the fields to the rectangular domain [—5, 5]^ and we 
consider a two-dimension uniform-space discretization, namely, Xkj = (—5 -I- k6,—5 + jS), where 
S = 10/(2m — 1), for k,j = 0,..., 2m — 1. We use m = 128. 

The values of the norms of the scattered electric and magnetic fields |mo|) I^^oI and the interior 
electric and magnetic fields |ui|, |ui| are presented in Figures 4 and 5 for different values of the 
polar angle </, which in corresponds to the incident direction (cos </, sin </) G S. 
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